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Abstract 

Results for the dynamic structure factor of a one-component plasma are presented for 
the situations of strong coupling and weak degeneracy and strong degeneracy and weak 
coupling. Possibilities to obtain rigorous results when both, coupling and degeneracy are 
large are discussed. 

In dense astrophysical and laboratory plasmas one frequently encounters situations 
where Coulomb and quantum effects are important simultaneously JT|, i.e. where the 
Coulomb coupling parameter T = (47rn e /3) 1 ' /3 e 2 /47reofcsT and the degeneracy param- 
eter x — n e^l [A e is the electron thermal wave length A 2 = 2nh 2 (3/m e \ exceed unity. 
While a rigorous theoretical treatment of the equilibrium properties is now becoming 
possible, in particular by path integral quantum Monte Carlo (PIMC) techniques [[|, 
comparable results for the dynamic properties are still missing. In this paper, we dis- 
cuss two possible approaches to solve this problem. We focus on the dynamic dielectric 
function e(u>, q) and the dynamic structure factor S(ui, q) ~ — Im l/e(u>, q) from which 
all optical and transport properties can be computed. 

The central quantity which describes dynamical properties of an N-particle system 
is the density fluctuation p q , where q is the wave vector. In a classical system, it can 

be computed from Klimontovich's phase space density p p>r (t) = £ i=1 $l r ~ r i(*)] ^[p ~~ 
Pi(£)], where Ti(t) and Pi(t) are classical (random) particle trajectories. Fourier trans- 
formation with respect to r yields the momentum- dependent density fluctuation p p . q , 
and the total fluctuation p q is obtained by a subsequent integration over all momenta, 

N N 

PpM = J>^) S[p - Mt)}, p q (t) = £ C *A<*>. (1) 

i=l i=l 

The corresponding quantum result is derived from the N-particle wave function 

if N {v,t) = £ a a a (*)0#(r,*), where is a complete set of basis functions, and 
the coefficients are creation/annihiliation operators obeing ( ant i-) commutation rules 
[al(t),d s i(t)] T = 5 S ^ S > for bosons (fermions). The quantum generalization of Eq. (|I|) is 
then given by (h = 1) 

Pp,qW = « P (*2) a P +q(*l)l tl =t 2 = t > Pq(*) = ^ ^(t) a p+q (t), (2) 
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where the fluctuation p P!q describes the annihilation of an electron in momentum state 
p + q at time t\ and creation in momentum state p at time t%. Obviously, the limiting 
case q = of Eqs. (JTJ) and (0) yields the momentum distribution and the density. 

The fluctuations p P;q and p q are random quantities depending on the initial condi- 
tions of the equations of motion. To come to measurable quantities, a suitable average 
over the particle ensemble (e.g. over the initial conditions) has to be performed yield- 
ing p p ,q(t) = (Pp,q(t)) and Pq(t) = (p q (t)). These are already macroscopic observables 
which are closely related to the dynamic structure factor. There exist two ways to 
compute S(u>,q) for a N-particle system in equilibrium. 

The first is to analyze the equilibrium fluctuations in a field-free system. In 
this case, p q (i) ~ due to spatial homogeneity. The first non- vanishing average is the 
density-density correlation function C pp which is closely related to pair correlations 
and the Fourier transform of which yields the dynamic structure factor, 

1 f°° 

C w (q,T) = (p q (r)p_ q (0)), fif(w,q) = — / dr C w (q, r). (3) 

Cpp can be computed from equilibrium theories, such as the Bethe Salpeter equation, 
or nonequilibrium approaches, in particular molecular dynamics simulations (MD). 
The latter has the advantage that, for a classical system, correlations of arbitrary 
strength are straightforwardly included, e.g. in the frame of Newton's equations. After 
averaging the resulting trajectories fj(t) and p,(i) [over runs with different initial 
conditions or over time], the dynamic structure factor follows directly from Eqs. ([I]) 
and ©. 

We have applied this approach to a degenerate electron gas by using an effective 
quantum pair potential in the classical MD scheme. An appropriate potential was ob- 
tained by Kelbg 0] and is rigorous for V < 1 and \ < !• We have performed extensive 
simulations in which the influence of quantum effects on the dynamic structure factor 
and the plasmon spectrum has been investigated 0. It is found that quantum ef- 
fects influence the plasmon spectrum predominantly at large wave numbers, but their 
influence is rather small. Fig. l.a shows an example of the simulation results and a 
comparison to mean-field theories. 

The second approach to dynamic properties is to investigate the temporal re- 
sponse of the plasma to a weak inhomogeneous field U(t) with wave number 
q Q . The field gives rise to collective oscillations of the electrons and thus to a non- 
vanishing fluctuation p qo (t) the Fourier transform of which directly yields the dynamic 
structure factor |3j 

S(«,q.)~-i^M. (4) 

This expression is much simpler than (|3]) and well suited for classical MD simula- 
tions. On the other hand, using the quantum expressions (||]), a kinetic theory of 
the dynamic properties of strongly degenerate systems can be developed. Starting 
from the Heisenberg equation of motion % ^p P , q — [p p , q , H sys + H ext ] = 0, with the sys- 
tem and external hamiltonians H sys = J2k + \ Z]kik 2 q^o ^ / (?)^k 1 + q "k 2 -q^ k 2^ki, 
H ext (t) = Eq^^i ^)Ek4^k-qi one can derive an equation of motion for the av- 
eraged fluctuation of p p , q . However, the most general theory is obtained if one 
considers the two-time fluctuations, cf. Eq. (0), [there are 2 independent functions 
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due the non-commutivity of the operators a,af] Gp liP2 (£i£2) = *(o P2 (t2)api(£i)) and 
G>„ (^1^2) = — *(o pl (ii)a Pa (i2))- The equations of motion for are the inhomoge- 
Kadanoff-Bavm eauations fKBE") and read 



Pl,P2 

neous 



v-j-/-p 2 --a 

Kadanoff-Baym equations (KBE) and read 
1 q 

E E S(*^)G ! f ik2 (^2)+/f 1)k2 (tlt2), 



(5) 
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and are to be supplemented with the adjoint equation. The Coulomb potential V 
gives rise to mean-field and short-range correlation effects which are contained in the 
Hartree-Fock selfenergy H HF and the collision integrals I, respectively. 

Eqs. (U) are very gen- 
eral and include the re- 
sponse from a nonequi- 
librium state as well as 
the nonlinear response in 
case of a strong external 
field. But here our inter- 
est is in the linear response 
to a weak field resulting 
in Eq. (f|). This limit 
is obtained by considering 
Eqs. (pj) for a monochro- 
matic excitation £/(q, t) = 
£/o(t)<5 q>qo of small ampli- 
tude Uo which allows for 
a linearization. The result 
are coupled equations for 
the spatially homogeneous 
and monochromatic com- 
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Fig. 1: Dynamic structure factor for an electron gas. a: 
Using Eq.(3) with molecular dynamics simulations with the 
Coulomb and Kelbg potential compared to the RPA and 
Vlasov theories for weak degeneracy, nA 3 = 0.5, T = 1, 
qf = 0.6. b: Quantum kinetic theory results for S, Eq.(4), 
comparing RPA, RPA plus exchange and RPA with exchange 
and collisions in Born approximation for strong degeneracy, 
r s = 4, T = O.IQEji, q = 0.6/cp. f,r s ,En and hp denote the 
mean interparticld distance, Brueckner parameter, Hydrogen 
binding energy and Fermi momentum, respectively. 
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l^iol 



Uq, for details, 



cf. Ref. 0. After solving this system, the total fluctuation follows according to 
Pq o C0 = — iJ2 p Gp +Ci ;P (t,t), the Fourier transform of which yields the structure fac- 
tor (£|). This is illustrated in Fig. 2, middle and lower figures. 

Equations (Eh contain all important limiting cases of a quantum electron plasma: 
i) neglect of S and I yields just the spectrum of single-particle excitations; ii) 
inclusion of T, H (Hartree mean field) is equivalent to the RPA (or - in the classical 
limit - the Vlasov response); iii) inclusion of S HF = + yields the RPA corrected 
by exchange (time-dependent Hartree-Fock) and, finally, iv) addition of I< yields RPA 
with exchange and correlations. Numerical results comparing cases ii)— iv) are shown 
in Fig.l.b. 

Most remarkably, consistency problems (e.g. sum rule preservation) often encoun- 
tered in theories incorporating collisions in the dielectric response functions, are com- 
pletely avoided due to the conservation properties of the Kadanoff-Baym equations. 
In fact, choice of density conserving collision integrals is sufficient to preserve the f- 
sum rule. Moreover, the KBE easily allow to fully conserve the total energy, i.e. all 
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00 , ^r 10 are full renormalized Greens functions and are thus appropriate 



5 0.75 
| 0.5 



T=0.16E R 



4 6 8 

wave number k, l/a B 



functions Gj§ 

to describe the response from a correlated equilibrium or nonequilibrium state of the 
plasma. Interestingly, the solution of the KBE in the time domain allows to account 
for very complex correlation effects by 
using very simple (compared e.g. to the 
Bethe Salpeter theory) selfenergies in 
the collision integrals, and it automati- 
cally yields all necessary vertex correc- 
tions. The numerical result in Fig.l.b 
used the static Born approximation, 
thereby fulfulling the f-sum rule to bet- 
ter than 0.1%. 

To extend these calculations to 
plasmas with strong correlations and 
strong degeneracy, one has two possibil- 
ities: 1. Generalize MD simulations to 
nonequilibrium quantum systems, e.g. 
[[J. 2. Include strong correlation effects 
into the collision integrals of the quan- 
tum kinetic equations. Alternatively, 
one can use, instead of the selfenergies, 
pair correlation functions from inde- 
pendent calculations M, such as path 
integral Monte Carlo . This should 
allow to obtain first-principle dynamic 
response functions in the near future. 
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Fig. 2: Quantum kinetic results including RPA 
and exchange for a strongly degenerate electron 
gas in an optically excited n-doped semicon- 
ductor, upper fig.: momentum distribution, 
with equilibrium (density n E ®, temperature T) 
and nonequilibrium electrons (n NE ). Middle 
fig.: time-dependent total density fluctuation 
(full line) excited by a monochromatic external 
field (dashed line). Lower fig: structure factor 
with an optical (high frequency) and a weakly 
damped acoustic plasmon which emerges out of 
the pair continuum. 
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